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Abstract 

We present exact inhomogeneous and anisotropic cosmological solutions of low-energy string theory 
containing dilaton and axion fields. The spacetime metric possesses cylindrical symmetry. The solutions 
describe ever-expanding universes with an initial curvature singularity and contain known homogeneous so¬ 
lutions as subcases. The asymptotic form of the solution near the initial singularity has a spatially-varying 
Kasner-like form. The inhomogeneous axion and dilaton fields are found to evolve quasi-homogeneously 
on scales larger than the particle horizon. When the inhomogeneities enter the horizon they oscillate as 
non-linear waves and the inhomogeneities attentuate. When the inhomogeneities are small they behave 
like small perturbations of homogeneous universes. The manifestation of duality and the asymptotic 
behaviour of the solutions are investigated. 

PACS numbers: 98.80Hw, 04.20.Jb,04.50.+h, ll.25.Mj 


1 Introduction 

The low-energy effective action of the bosonic sector of string theory provides a gravitation theory containing 
dilaton and axion fields that possesses cosmological solutions. These solutions provide models for the 
behavior of the universe near the Planck (or string) energy scale [1]. They allow us to investigate a number 
of long-standing questions regarding the occurrence of singularities, the behavior of the general solution 
of the theory in the vicinity of a singularity, and the likelihood of our Universe arising from generic initial 
data. They also provide a basis for investigation of higher-order corrections to low-energy cosmological string 
theory. Several studies have recently been made of string cosmologies in order to ascertain the behavior of 
simple isotropic and anisotropic universes, investigate the implications of duality, and search for inflationary 
solutions [2]-[8]. Since one of the prime reasons for studying such solutions is to shed light on the behavior of 
the universe at very high energies, where our knowledge of its material content, geometrical and topological 
properties, or its anisotropies and inhomogeneities, is necessarily incomplete, it is unwise to make special 
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assumptions about the form of the cosmological solutions. Indeed, any dimensional reduction process could 
be viewed as an extreme form of anisotropic evolution in more than three dimensions in which three spatial 
dimensions expand whilst the rest remain static. A number of studies have focused on obtaining particular 
solutions for 3+1 dimensional spacetimes in cases where spatial homogeneity (and sometimes also isotropy) 
is assumed for the metric of spacetime, where the H field is set to zero [4], or where the H field is included 
by assuming that it takes a particular form which satisfies its constraints and its equation of motion [5]. 
For example, Copeland et al, [2], discussed Friedmann and Bianchi type I universes, allowing *H to be 
time dependent or space dependent, respectively. In a second paper, [3], they discussed Bianchi I solutions 
with a homogeneous antisymmetric tensor held. In [6] (see also [5]) Batakis presented an overview of all 
possible configurations of a (spatially) homogeneous IF-held in diagonal Bianchi models. Whereas, in ref. 
[7], we investigated the case for a (spatially) homogeneous tensor potential B^ in Bianchi metrics that 
are not necessarily diagonal. We also gave a classification of all the degrees of freedom permitted for the 
H held in spatially homogeneous universes possessing a 3-parameter group of motions. The only spatially 
homogeneous universe excluded from this study is the (closed) S 2 x S 1 Kantowski -Sachs universe. A detailed 
study of this universe was made by Barrow and Dabrowski [8]. 

In this paper, we take one further step upwards in generality and consider a wide class of inhomogeneous 
and anisotropic string cosmologies. These possess cylindrical symmetry and contain homogeneous Bianchi 
and Kantowski-Sachs universes as special cases [9]. They allow us to investigate the propagation of non-linear 
inhonrogeneities in the axion and dilaton helds. On scales larger than the particle horizon inhonrogeneities 
in the axion and dilaton helds evolve quasi-homogeneously but when the inhomogeneities enter the horizon 
they undergo oscillations and attentuate. In the limit that the amplitude of the inhonrogeneities is small we 
will recover the results of perturbation studies of homogeneous string cosmologies in an appropriate gauge. 
Besides providing exact descriptions of the gravitational self-interaction of strongly inhomogeneous axion 
and dilaton helds these solutions allow us to investigate the impact of duality upon the form of the solution 
in a situation where there exist characteristic spatial scales. 

The string world-sheet action for a closed bosonic string in a background held including all the massless 
states of the string as part of the background is given by, [1], 

5 = Vhh a ^d a X p d p X v g^{X p ) + e a0 d o X p dgX u B^(X p ) + a'Vh<t>{X p )R {2) } (1) 

where h a @ is the 2-dimensional worldsheet metric, R 1 ' 2 ’ 1 the worldsheet Ricci scalar, e“ /5 the worldsheet 
antisymmetric tensor, B^ IU (X P ) the antisymmetric tensor held, g^ lu (X p ) the background spacetime metric 
(graviton), 4>(X P ) the dilaton, a! is the inverse string tension, and the functions X p (a) map the string 
worldsheet into the physical D-dimensional spacetime manifold. 

For the consistency of string theory it is essential that local scale invariance holds. Imposing this 
condition results in equations of motion for the helds g /w , £> /i;/ and 4> which can be derived to lowest order 
in a! from the low-energy effective action for a vanishing cosmological constant 

S = J d D x^ge~\R + g ap d^dpcf) - ^H a ^H a01 ). (2) 
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where the antisymmetric tensor field strength = d\ a B^ 1 is introduced. 

In a cosmological context it is generally assumed that by some means all but four of the 10 or 26 
dimensions of spacetime are compactified, leaving an expanding 3+1-dimensional spacetime (D = 4). Since 
we are interested in cosmological solutions of the field equations derived from the variation of this action, 
we adopt the Einstein frame by making the conformal transformation 

9a0 * C 9a0- (3) 

In this frame the 4-dimensional string field equations and the equations of motion are given by (Greek 
indices run 0 < a, (3 < 3). 

The low-energy effective action in the Einstein frame yields to the following set of equations {k 2 = 
8irG, c = 1) 

R^-\g^R = k^t^+Wt^) 

= 0 

□< 3 ? + = 0 

where 

(7) 

(H % V = ^e- 2$ (3 (8) 

Thus in this frame the problem reduces to the study of inhomogeneous general relativistic cosmologies 
containing two gravitationally interacting matter fields. In the next section we shall manipulate these 
equations into a soluble system by introducing a particular inhomogeneous spacetime metric with cylindrical 
symmetry with a particular topology. In section 3 we give exact solutions in cases where one (or both) of 
the axion and dilaton fields depends only upon the time variable. In section 4 we consider the case where 
both fields depend upon time and space coordinates. In section 5 we investigate the asymptotic behaviours 
of these fully inhomogeneous solutions on scales large and smaller than the horizon. In section 6 we study 
the nature of duality in these solutions and the results are discussed in section 7. 

2 Einstein-Rosen Metric 

Consider the anisotropic and inhomogeneous spacetime metric [10] 

ds 2 = - e 2 (x-^) ( dr 2 - dr 2 ) + R(e 2 ^dz 2 + e" 2 ^# 2 ) (9) 


(4) 

(5) 

( 6 ) 
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where Xi'&iR are unknown functions of t and r. Thus d/dz and d/d<f> are Killing vectors. Without loss of 
generality it can be assumed that 0 < z < 1 and 0 < (f> < 1. When ij) = 0 and R = e 2x , with R = R(t) 
and x = x(t), we recover an isotropic Friedmann universe. Other homogeneous specialisations of the metric 
reduce it to one of the Bianchi type homogeneous universes [ 9 ]. Properties of the metric ( 9 ) depend on 
whether B^ = R ^ is spacelike, timelike or null (Greek indices run 0 —> 3). The cases with a globally null 
or spacelike B^ correspond to plane or cylindrical gravitational waves, respectively [2], Metrics where the 
sign of B^B 11 varies throughout the spacetime describe colliding gravitational waves [3] or cosmologies with 
timelike and spacelike singularities [4], Metrics with a globally timelike B M describe cosmological models 
with spacelike singularities. If the spacelike hypersurfaces are compact then the allowed spatial topologies 
[5], [ 9 ] are a 3-torus, S 1 <g> S ' 1 <8> •S' 1 , for R = (detg a b) 2 = r with 0 < t < 00 and 0 < r < 00; a hypertorus, 
S ' 1 <g> S' 2 , or a 3-sphere, S 3 , for R = (detg a b) 2 = sin r sin r with 0 < r < ir and 0 < r < ir. We shall present 
solutions for the globally timelike case, R = r. These correspond to ever-expanding cosmological models 
with an initial curvature singularity at r = 0. Note that the behavior of the closed S 3 models approaches 
that of the 3-torus universes as the singularities are approached because sinr —* r as r —> 0 and 7 r, and so 
the role played by the duality invariance of these models can be investigated along with the implications for 
the ’pre big bang’ scenario of Gasperini et al [4]. The homogeneous models of the S 1 <g> S 2 case will be the 
Kantowski-Sachs universes studied in ref. [8]. For further results about the singularity structure and global 
existence of these metrics (the strong cosmic censorship hypothesis holds) see the paper by Chrusciel et al 

[15]- 

Rewriting (5) as 


d(*H) - 2(d$) A (*H) = 0, 


( 10 ) 


and using 


dH = 0, 


( 11 ) 


we can determine the general form of H that is compatible with the Einstein-Rosen spacetime geometry. 
Denoting x° = r, x 1 = r, x 2 = z, and x 3 = 4> we require 

H = 6^4(t, r)dx° A dx 1 A dx 2 + QB(t, r)dx° A dx 1 A dx 3 (12) 

+6C(r, r)dx° A dx 2 A dx 3 + 6 D(t, r)dx x A dx 2 A dx 3 . 

The quantities H and can be a functions only of r and t here since the energy-momentum tensor is 
allowed to depend only on these variables. Hence (11) implies 

d 1 C-d 0 D = 0, (13) 


while -kH is given by 


( 14 ) 


4 



Since d<& = do&dx° + di&dx 1 , eqn. (10) reduces to 

[<9 0 Fi - 3iF 0 - 2(Fid 0 $ - Fod^jdx 0 A dx 1 
+ [d 0 F 2 - 2F 2 do<f>]dx° A dx 2 + [d\F 2 - 2F 2 d l <S>]dx l A dx 2 


+ [doF 3 - 

2Fodo&]dx° A dx 3 + [<9iF 3 — 2F 3 (9i < 3?]dx 1 A dx 3 = 

0. 

(15) 

This implies 

doF\ — 2F\do < & — d±Fo + 2Fo<9i < h = 0 


(16) 


d{e~ 2q, F 2 ) = 0 


(17) 


d(e~ 2 *F 3 ) = 0, 


(18) 

so (17) and (18) yield 

F 2 = e 2 ®A 2 , 


(19) 


F 3 = e 2 ®A 3 , 


(20) 

where A 2 and A 3 are constants 




Using the fact that C(r,r) 

= 9m922g 33 H Q23 = goo922933£ 0231 Fi = —RF\ and, 

similarly that D(r, 

r) = 

— RFo, eqn. (13) becomes 

di(RF) - do(RFo) = Q, 


(21) 

and (16) implies that 

doie-^h) = d 1 (e~ 2 *Fo). 


(22) 


In order to solve the system of differential equations (21)-(22) there are two obvious choices: 

(i) e^'^Fi = d\b and e~ 2iS> Fo = dob 

(ii) RF\ = doh and RFq = d\ h 


The latter choice corresponds to taking B 23 to be the only non-vanishing component of the antisymmetric 
tensor potential defined by H = dB and depending only on r and r. 

The choice (i) reduces (21) to 


□6 + 2V Al 6V Al <h = 0, 

while choice (ii) produces another coupled wave equation 

h - h" ~^h+ - 2(<M - $'/;,') 

IX IX 


= o 


(23) 


(24) 
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where • = and ' = §p. 


The off-diagonal components of the Einstein tensor G 02 , G 03 , G\ 2 , G 13 and G 23 , are zero in the spacetime 
(9). The corresponding components of all vanish so we only need ensure that all the corresponding 

off-diagonal components of are also zero. Since we have 

(H) T 02 = - 2 F 2 F 0 e~ 2 * = l -A 2 F Q , (25) 


we must therefore set A 2 = 0, and hence H 013 = B = 0. For the (03) component we have 



{H) T m = ^Tbe" 2 * = ±A 3 F 0 , 

(26) 

so we must set A% = 0, and hence H 120 = A = 0. With these choices, the components t ' H )T\ 2 , and 

(^)?23 also all vanish. Therefore, the equations governing the dilaton and antisymmetric tensor field for the 
two choices are given by equation ( 6 ) together with the following coupled propagation equations, (27)-(28) 
or (29)-(30), in the cases (i) and (ii), respectively: 

' (i) 


4 (#$')'-4 (it®) ~e 2 *\b' 2 -b 2 1 = 0 , 

(27) 


(Rb')' - i (m>y + 2 [$v - <±>4 = 0 . 

(28) 

(ii) 

R-\m'y - R-\R®y + R- 2 e- 2q, [ti 2 - h 2 } = 0, 

h - h" - 4 h + ^-ti = -2 {Q'ti - 6 / 1 ). 

R R 

(29) 

(30) 


Since both choices involve the same number of independent functions they are equivalent; here, choice 
(i) is taken. 

The energy-momentum tensor in (4) reads 


k 2 (*) t ^ = l _yx d ^ d ^ _ I«S;)(3<1>) 2 ), 

(31) 

(H)y = \y {g ux d ^ bdvb _^_ 5 x {db) 2 l 

(32) 
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so the nonvanishing components of the energy-momentum tensor are 


= [T 2 + 2 + {b 2 + b' 2 )e 2<1> ] = -k 2 ^ 1 , 

= i e - 2 (x-i« [$$' + M/e 2 *] = -k 2 ^ 0 , 

= i e - 2 (*-^) [<F 2 - $' 2 + (fe 2 - 6' 2 )e 2<J> ] = k 2 T 3 3 , 

(33) 

The energy-momentum tensor can be interpreted as describing two stiff perfect fluids where the energy 
density for the dilaton fluid is found to be 


« 2 V 

k 2 T 2 2 


p$ = / ^$ = ie- 2(x - 1/ ' ) [i 2 -^> ,2 ] 

and its 4-velocity is given by 

u a = e -(x—0) [<j,2 _ 2j-| 3/ 0 , o) 

and for the axion fluid we have 

Ph = PH = ^e~ 2 ^ x ~^[b 2 — b' 2 ]e 2$ 

and its 4-velocity is 

v° = e~^-^ [b 2 - b' 2 }~H-b, 6', 0,0). 

Furthermore, T 0 0 + 1 = 0 and T 2 2 — T 3 3 = 0 as the cylindrical symmetry of the metric demands, and 


Einstein’s equations for i?(r, r) and ip(r, r) are given by the linear wave equations [11] 

R-R" = 0, (34) 

$ + ^ ~ ~ = 0. (35) 

The remaining metric function, y(r, r), is determined by the two Einstein constraint equations 

x' = $ - (R 2 -R ,2 r 1 [RR\'>P 2 + ^' 2 ) + R'R"-ZRR^'-RR' (36) 

4 ix 

-K 2 Re 2{x ~^ (T 0 °R ' + T 0 1 R)], 

X = ip-~-(R 2 -R' 2 )- 1 [2RR , iP'ip , -RR('ip 2 + il> ,2 )-RR" + R , R' (37) 

4 ix 

+n 2 e 2 ^- i ’ ) R{T 0 °R + T 0 ^O]- 
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Since cosmological solutions are of interest to us here, we consider only the timelike solution of (34). 
Using the general coordinate invariances (r ± r —► f(r ± r)) of the metric this may be taken without loss of 
generality to be 


R = R(t) = t. 

Then, equation (35) reduces to 

Pp + —ip — ip" = 0 

T 

which is solved by 


(38) 


(39) 


ip(T,r) =° ip 0 +° V’l lnr + cos K r - r n )][ A ^ n J 0 {riT) + B T„iVo(nr)] (40) 

n =1 

where °ipi, B ^ n , r n are constants and Jq(x) and TVq(x') denote the 0 th order Bessel and Neumann 

functions, respectively. 

Equations (27)-(28) read 


-$- e 


2$ 


b' 2 — b 2 


b” - 6 + 2 


$' 6 ' - $6 


In the next section several solutions will be found. 


(41) 

(42) 


3 Solutions of varying generality 


Before explicit solutions are given something general we make some remarks about the procedure for solving 
the system of partial differential equations for the metric function y(r, r). 

For R(t) = r, eqns. (36) & (37) reduce to 

X 1 = ip' + 2npip' + K 2 re 2 ^ x ~^T 0 1 (43) 

x = V’ —-— + Tpip 2 + ip' 2 ) — K 2 re 2 ^ x_ ^T 0 0 (44) 

4 T 

Generally speaking, the most difficult step is to find the integral for the part coupled to ip(r, r ). However, 
this problem was solved by Charach [16]. Define a function G(ip;r,r) by 


G' = ripip' 

= ^r(ip 2 + ip' 2 ) 
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G 


(45) 

(46) 



Note that ip satisfying 


-0 H —ip — ip" = 0 
r 

is kept as a functional dependence in G. The explicit dependence on t and r might sometimes be suppressed 
and we write G(ip) as 

G(ip',r, r) = Vo + 7j (W^ir 

OO 

+Vi 51 cos H r - ’’n)][ A ^ r n^o(^) + B ^ n N 0 (nr)] 

n=l 
■ oo 

+-T 2 55 n 2 ([' 4 T n Jo(nr) + B ^„AT 0 (nr)] 2 + [ j4 T„Ji(nr) + B ^ n 7Vi(rar)] 2 ) 

4 n=l 

^ OO 

--r 55 ncos 2 [n(r - r n )]{( 4 4/ n ) 2 Jo(nr) Ji(nr) 

2 n=l 

+ A ^^ n [N 0 (nT)J 1 (nT) + J 0 (nr) A^i (nr)] + ( i3 'I' n ) 2 iVo(nr)A r i(?rr)} 

-j OO OO 

+ 9 r 5Z 5Z 2 { sin [ ra ( r “ V] sin[m(r - r m )][nC/^,(^) - ™€( T )1 

Z mi/ 

n= 1 m=l,n0TO 

+ cos[n(r - r n )] cos[m(r - r m )][m[/f2,(r) - nt4m( r )]} 

(47) 

where 


U ml( T ) = A V )Jo(nT ) + B ^ B ^ m iVo(rnr)iVi(nr) + 2 j4 T B T m Ji(rar)iV 0 (mr), 

V^(r) ^ 4 ^ m Jo(nr)Jo(mr) + B ^ B ^ m iV 0 (nr)JVi(mr) + 2 y V B ’I' m J 0 (nr)JVi(mr). 

We now consider classes of solutions in which one (or both) of the <J> and b fields depend on only one of 
the coordinates r and t. 


3.1 Solutions homogeneous in r : $ = <3 >(t), b = b{r ) 
The well-known solution to (41) & (42) [6] in this subcase is 


e 


<i> 


K 0 


/ B 2 

cosh(iVC) + y 1 ~ ^2 sin h(iVC) 


jysinh(jVC) + ^/l-||cosh(jVC) 
5 cosh(lVC) + ^/l - sinh(JVC) 


(48) 

(49) 
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where N and B are constants and dr = rd(. 

Using this in the expression for the components of the energy-momentum tensor gives an expression for 
X(r,r) 


x{t, r) = ip(r, r) + 2G(ip; r, r) 


N 2 - 1 

---In r + M, 

A ‘ 


(50) 


where M is a constant. Hence, the metric function exp[y — ip] is given by 


exp[y -ip]= e 2G ^)e M r4 i (51) 

3.2 Solutions homogeneous in r : $ = $(r), b — b(r) 

The solution in this subcase is given by 

I B 2 

e® = cosh(Nr) + y 1 — sinh(Nr) (52) 

N sinh (Nr) + cosh (Nr) 

b(r) = —- = - (53) 

cosh(IW) + yi - ^ sinh(Nr) 

where N, B constants. From these expressions, y(r, r ) is found to be 


1 N 2 

X(r, r) = ip(r, r ) + 2G(ip; r, r) - - In r + — r 2 + M, 
which gives the remaining metric component, 


(54) 


exp[y - ip] 


= e 2G^) e 2M T - 


1 N 2 _2 

4e s 


(55) 


3.3 Solutions with an oscillatory axion : <4> = $(r), b = b(r,r) 

If we rewrite (42) as 

d 2 b d 2 b db d§ 1 ^d& db 
dr 2 8 t 2 dr "" dr r + ^ dr dr 

and take a solution 24>(r) = — ln(r/ro), then the axion field b(r, t) also satisfies the wave equation, 


(56) 


d 2 b d 2 b 
dr 2 dr 2 


(57) 
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which has the general solution 


(58) 


b(r, t) = abi(r + r) + /36 2 (r - r) 

where a, [3 are constants and are arbitrary functions of their arguments. Eqn. (41) is satisfied if 



This implies 

/ db db\ (db db\ 

V dr dr ) \ dr dr ) 

so that either a or (3 must vanish. Thus we obtain the solution 

T(r) = — ^ In —, 

1 T 0 

b(r,r) = Q(u)bi(r + r) + (1 - @(«))6 2 (t - r). 


(59) 


(60) 


(61) 

(62) 


where 0(it) is the step function (0(tt) = 0 for u < 0; 0(it) = 1 for u > 0) and u an arbitrary real parameter. 

It is interesting to have a solution with a homogeneous dilaton and an inhomogeneous axion. Note, that 
in this case the axion behaves quite differently from the dilaton. 

The function y(r, r) is determined by 


X 1 = V’ / + 2 Tipij:' + y bb' 

X = ^--|-+T(^ 2 +^ 2 ) + ^(fe 2 + 6 /2 ). 

16r 4 

To solve this system of equations we define a new function B(r,r) satisfying 


Changing to new variables, 


we find 


This implies 


B' = bb' 

B = l -{tf + b' 2 ) 


X = r + t , Y = r — t, 


dB 

f db \ 2 

dB 

fdb\ 

dx ~ 

laiJ ’ 

dY ~ ~ 

\dY ) 


d 2 B 

= 0 



dXdY 



(63) 

(64) 


(65) 

( 66 ) 


( 67 ) 
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which is generally solved by 


B(X,Y) =B l {X) + B 2 (Y) (68) 

with Bi arbitrary functions of their arguments. Using the general solution for b in terms of X and Y, 
B(X,Y) is given by 


B(X,Y) = @(u) j dX - (1 - ©(«)) J dY (f) 2 (69) 

Finally, an expression for %(r, r) is obtained, 

X(r, r) = r) + 2 + |fl(r, r) - ^ Inr + M, (70) 

which results in 

exp[x -i/j\ = e 2G W e M T -T6 e Wr,r)_ (71) 

3.4 Solutions with = 4>(r, r), 6 = b(r) 

If we take b(r) = At 2 /2, A constant, and 4>(r, r) = — lnr + S(r), then eqn. (41) requires S(r) to satisfy 

Hence, 

e -$(r,r) 

b(r ) 

is a solution of (41)-(42). 

Calculating the appropriate components of the energy-momentum tensor yields 


dr 2 


+ A 2 e 2S ^ = 0. 


/ A 2 

r[cosh(lVr) + y 1 — sinh(lVr)], 
^r 2 . 


(72) 

(73) 


x(r, r) = V’(u r) + 2r, r) + - ln[cosh(lVr) + 


/ H 2 j\r 2 

1 - ^ sinh(fVr)] + —r 2 + M, 


hence 


(74) 


exp[x — V’] = e 2G ^e M [cosh(iVr) + 


/ H 2 

1 — -^2 sinh(.ZW)] 


1 JV 2 - 2 
2e 8 


(75) 
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3.5 Discussion 

Apart from case (3.3) the solutions presented so far describe non-oscillatory axion-dilaton systems on an 
oscillatory cosmological background. In case (3.3) the axion field is allowed to oscillate which couples the 
dilatonic and gravitational waves. However, because of condition (59), only travelling wave solutions in 
6 (r, r) are described in this case. 


4 Charach Solutions 


The system of equations (41)-(42) is very similar to equations determining the components of the electromag¬ 
netic potential in the electromagnetic Gowdy universe [16] [17]. It was stated in [16] (and references therein) 
that the geometric requirements of the Einstein-Rosen spacetimes allow four independent components of the 
6 possible components of the Maxwell tensor which can be derived from two non-vanishing components of 
the electromagnetic potential. In section 2, we found that only two of the four possible components of the 
antisymmetric tensor field strength can be non-vanishing, which can then be accordingly derived from the 
potential like function 5(r, r) or /i(r, r). In the latter case there is a direct connection to the antisymmetric 
tensor field potential B^ u where H = dB. In order to obtain an exact solution of (41)-(42) where $ and b 
are dependent on r and r we employ a procedure introduced by Charach [16]. 

Assume that 

4>(r, r) = — - In v[b(r, r)] (76) 


where v(b) is a function yet to be determined. Since 

• 1 ■ d In v 

4> = --b ——, 

2 db ’ 

- 1 --d In v ■ n d 2 In v , 


= _l,dlnv 
2 db ’ 

I un d hi v . 2 d 2 In v 

* =-2 [b \ar + b 


we can use (42) to transform eqn. (41) into 


'd 2 v 


+ 2^ (V 2 - 6 2 ) 


= 0 , 


while eqn. (42) becomes 


6"-s 

r db 


(77) 


(78) 


Assuming b 1 “ — b 2 / 0, eqn. (77) implies 


v(b) = -b 2 + ci b + c 2 , 


(79) 
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where the logarithm in (76) requires that the constants ci, C 2 satisfy the inequality 


cf + 4c 2 > 0. 


( 80 ) 


Equation (78) becomes 


b" -b- 



T 


Cl , 2 (b 1 2 — b 2 \ 
c 2 + ci b — b 2 V / 


(81) 


If we make the substitution 


b = b{) + M tanh(Mcj) 


(82) 


eqn. (81) becomes 
M 2 


cosh “(My 


■ [u" — Co - Co\ + 2 


M 3 siiih(Mct-’) r . 
cosh 3 (Mw) 


- w' 2 1 = 2 


) Af 3 sinh(Aftu) f . 
cosh 3 (Aiw) 


tc 2 — U / 2 1 


where 6 q = 4ci and Af 2 = c 2 + 4c 2 > 0. Hence, co(r, r) satisfies the linear wave equation, 


(83) 


U> — LO - 00 = 0. 

T 


The wave-packet solution to (84) is given by 

OO 

uj(t, r) =° LOO +° CJ 1 lnr + ^ cos[n(r - r n )][ A H n J 0 (nr) + B H n lY 0 (nr)] 


71=1 


where °tUj, and r n are constants. 

In summary, equations (41) and (42) admit the following inhomogeneous solution: 


4>(r, r) = In 

1 


cosh(Afw) 

M 


b(r,r) = -ci + M tanh(Aftc). 

Rewriting the components T 0 ° and T 1 ° of the energy-momentum tensor in terms of co(r,r ) gives 

K VH) ro ° = _^ 2 + w ' 2 ), 


(84) 


(85) 


( 86 ) 

(87) 

( 88 ) 
(89) 
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and equations (43) and (44), which determine x(r, r), reduce to 


x' = V’ 7 + 2 Tipi/j' H— —tlolo' 

X = ^ ~+■(/;' 2 ) + +oj' 2 ) 


(90) 

(91) 


Using the function G(f; r, r) where / + r 1 / — f" 
(90)-(91) lead to 


0 and G( f: r, r) is given by eqn. (47) and eqns. 


1 M 2 

dx = di/) — -d In r + 2 dG(i/}\ r, r) H— —dG(u;-, r, r) (92) 

which yields 

1 M 2 

X(r, r) = ^(r, r) - - In r + 2G(^; r, r) + — G(w; r, r) + L (93) 

where L is some constant. 

So that the metric function exp[y — is given by 

exp [x -i/>]= e L T“3 e 2G W +i ^ G H (94) 

5 Asymptotic behavior 

The existence of inhomogeneity in the solutions found in section 4 introduces characteristic length scales 
and the gravitational self-interaction of the dilatonic and axionic waves will differ over scales according as 
they are causally coherent or not. The horizon distance in the r direction is defined by ds 2 | z ,</>= 0; hence 

A r = f dr = t. (95) 

Jo 

Therefore, the combination nr in the solutions above can be interpreted as the ratio of the radial horizon 
distance to the coordinate wavelength, A, since n oc 1/A(n). There are two limiting cases to be considered: 
the case nr <C 1, when the comoving wavelength is much larger than the radial horizon scale, and the case 
nr> 1, when the wavelength of the inhonrogeneities is well within the horizon scale. We consider these two 
cases separately. 

The Charach solutions discussed in the last section are the most general ones of those given. Apart 
from solutions of section (3.3) the limiting properties of the other solutions are included in those of the 
Charach-type solutions. Therefore in this section only the asymptotes of these solutions, of section 4, are 
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discussed. Explicit formulae for the functions involved are given in the appendix along with some useful 
definitions. It is convenient to define metric functions 


Ai(t, r ) 

= exp[x(r,r) - ip(r,r)] 

(96) 

A 2 {t, r ) 

= T \ eV’lrD 

(97) 

A 3 (r,r) 

= T ^e~^ r,r ' > 

(98) 


5.1 The limit nr <C 1 

In this case 

Ai(r,r) ~ e L r _ 4 +272 ^ ;r ^ + ~ 72 ^ ; ^ exp[ 27 i(' 0 ; r) + ^- 7 i(ca; r)] (99) 

A 2 (T,r) ~ e“i ( 100 ) 

A 3 (r,r) ~ e -ai(^;r) T i -a 2 (^;r) ( 101 ) 

This limit corresponds to the case where the comoving wavelength is much larger than the (radial) 
horizon size or in other words the universe consists of causally disconnected regions. In this case one would 
not expect to have any oscillatory behavior in nr. 

Concentrating on the homogeneous limit for r approaching zero allows to discuss cosmological solutions 
near the singularity. The metric functions are found to approach 


Ai(t) 

„ t -H(Vi) 2 +T-(°«i) 2 

(102) 

4j(t) 

~ T 5+°* 

(103) 

A 3 (t) 


(104) 


Changing to proper time using, in the homogeneous limit, the relation 

t = J drA i(t), 

r(t) is found to be 

r oc il+(°v-i) 2 +4 i ( 0 “i) 2 . (105) 

Defining the Kasner-exponents p t , i= 1, 2, 3, by 

~ diag(- 1, t 2pi , t 2p2 , t 2p3 ) 
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they are found to be slowly spatially varying: 


Pi 


P2 


P 3 

They satisfy the algebraic constraints, 


4 + (%) 2 + ^(V) 2 
f + (°Vh) 2 + ^(W ’ 
i+Vi 

f + (Vi) 2 + ^(V) 2 ’ 
I- Vi 

f + (Vi) 2 + ^(V) 2 ' 



Erf 


i, 

1 M 2 (V ) 2 

2 [f + (Vi) 2 + ^(W] 2 ' 


( 106 ) 

(107) 

(108) 


(109) 

( 110 ) 


The fact that J2i=iPi < T where the equality holds in the vacuum case (M = 0), shows immediately 
that there are isotropic solutions. This feature is present in the matter-filled Gowdy solutions [16] [17] and 
in the spatially homogeneous Kasner universes containing stiff fluid. 

The axion-dilaton system is independent of the gravitational background in the sense that its determining 
equations (cf (41) and (42)) do not involve any of the metric functions apart from i?(r, r). However, due 
to the general structure of the equations the solutions for and b are very similar to those of the metric 
functions. As r —> 0, the dilaton and axion fields approach 


<h(r, r) 
6(r,r) 




r\j 


ln[ e Mai ( w;r ) T Ma2 ( a,;r ) -|- e ~Mcn(u>;r) T -Ma 2 (u>;rh _ } n 2M, 

^ ^2Mai (w;r)^.2Mo!2(w;r) _ ^ 

2° lJr M g2Moii (w;r)j2MQ2(w;r) i 


(in) 

( 112 ) 


We note that the early-time behaviour of these solutions falls under the category of ’velocity-dominated’ 
solutions used in studies of general relativistic cosmologies [15]. As the singularity is approached the spatial 
gradients become negligible with respect to the time derivatives, 3-curvature anisotropies are ignored, and 
velocities are assumed to less than the speed of light. This approximation does not encompass the most 
general known behaviour in general relativity, with the metric undergoing chaotic oscillations on approach 
to the singularity [24]. Chaos in string cosmologies will the subject of a separate study [25]. 


17 



5.2 The limit nr 3> 1 


In this case the comoving wavelength is smaller than the (radial) horizon size allowing interaction between 
different modes and hence an oscillatory behavior of the metric components. From the limits of i/’(r, r) and 
G(ip', t, r) given in the appendix it can be seen that iJj displays an oscillatory behavior whilst the oscillations 
in G are damped out. 


Ai(r,r) ~ e L r 4 + 2 '>' 4 W’)+ 2 74 ( a ') ex p[27 3 (y>) + ^-73(0;)] 

M 2 

x exp[(27 5 (^) + —75(w))r] (113) 

A 2 {r,r) ~ exp[T~^h(i/}] r, r)] (114) 

A^(r,r) ~ T 2~4 2 W’) e -/3 i W’) exp[—T _ 2/t(^;r, r)] (115) 

As can be easily seen from the definition of h(^; r, r) given in the appendix it satisfies the wave equation 
(in Minkowski space), 

h-h" = 0 . (116) 

The exponential in A\ ensures that the homogeneous limit is approached at large r, and is an anisotropic 
universe which can be at most axisymmetric (fa = 0). Since g^ v ~ diag(— A 2 , A 2 , A 2 , A 2 ) and, for 

large values of r, we have 

exp[2r ~2 h(tp; r, r)\ ~ 1 + r, r], 

and so g^ v can be written as the sum of a background part and a “wave” part h liu 

fljtv — h///' 4“ h/jiv: 

which are found to be 


^ = diag(-A 2 , A 2 , T '+W) e WiW y~W) (117) 

hfj, u = diag(0, 0, 2t^ + 2 ^ 2 ^ 1 ' > e 2 ^ 1 ^h{^\ r, r), — 2r2 _2 ^ 2 ^e~ 2 ^ 1 ^ ;t,T 4/i(^; r, r)) (118) 

The dilaton-axion system displays an oscillatory behavior as well, although as emphasized before, there 
is no interaction between gravitational and axion-dilaton waves. 

The asymptotes are given by 
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$(t, r) 


In 


e M/3i(^) r M/3 2 (u)) 


e -M/3i(aj) r -M/3 2 (^) 


+M [ e M M^T M ^)-h - e -«AH T -M&H4]| ! ( t(J ; Tir ) ] _ l n (2 M) 


1 e 2 Mft(ai) T 2Mft(w)[| ^ 2 J^f T zh(uj-,T, r)l — 1 

6(r, r) ~ -ci+M- T - 

2 e 2M/3i(u>) r 2M/3 2 M[i _)_ 2Mr~2h(u;T,r)} + 1 


(119) 

( 120 ) 


6 Duality 


By means of dimensional reduction we can show that the low-energy effective action (2) is invariant under 
global 0{d,d) transformations, where d < D refers to the number of coordinates it does not depend on 
[20]. So, if one assumes a spacetime of the form N x K where N is a {D — d)-dimensional spacetime with 
coordinates (y = 0,1, ..D — d— 1), and K a d-dimensional compact space with coordinates y a (a = 1, d), 
and furthermore that all fields are assumed to be independent of the y-coordinates of the “internal” space 
K, then using the notation of [20] we can rewrite (2) as 

5 = J^dx J K dyyf^e~\m) +rd^d^ - ( 121 ) 

The hatted quantities now refer to the .D-dimensional spacetime. Using the vielbein formalism, gpo is written 
as 


. t + A^Aty A ^ \ 
' l G at) ) 

where g is the metric on N and G a p the metric on K. 

Define a shifted dilaton 

$ = </>- i log detG Q p, 


( 122 ) 


(123) 


and a 2d x 2d matrix Q, written in d x d blocks, 


/ G^ 1 -G~ l B \ 

^ BG~ l G - BG~ l B J ' 


It can be shown that (121) is invariant under global 0(d,d ) transformations 


<1 6 Q -a QQQ t 


(124) 


(125) 


19 



where f7 £ 0(d, d), that is Q T r]Q = rj where 


Id 

0 


n = 


o 

id 


and I ( | is the d-dimensional unity matrix. 

In the case of a diagonal g^o and a vanishing R-field, with the choices d = D — 1 and 17 = r/, the scale- 
factor duality is recovered. This was first discussed by Veneziano [21]. In this case the duality transformation 
results in an inversion of the scale-factors in the string frame. For a comprehensive discussion of target-space 
duality see ref. [22], 

In the case of the Einstein-Rosen metric, (9), considered here the low-energy effective action (2) is 
invariant under 0(2, 2) transformations for D = 4. In this section the antisymmetric tensor potential B^ w 
is assumed to be vanishing. Transforming the metric (9) to the string frame, the “internal” metric G a g is 
found to be 


_ ( Re^e® 0 \ 

Ga/3 - { 0 Re-^e* ) 


(126) 


The shifted dilaton defined above remains invariant under 0(d, d) transformations and this implies that 
the dilaton itself transforms as 


, . 1, detG a s 

—► <1>-In- -v^-r 

2 detG d ^ al 

where the index “dtta/” indicates an 0(d, d) transformed quantity. 


(127) 


6.1 Generalized Scale-Factor Duality 


As mentioned above, picking = 7] results in the case of a diagonal metric depending on just one (time-like) 
coordinate and results in scale-factor duality. So choosing = ij in the Einstein-Rosen case, where there is 
dependence on time and space variables, could be called a generalized scale-factor duality. 

Transforming Q according to (125) results in 


so that 


which implies 


G-^G" 1 , (128) 

R dual e = irV^e-* (129) 

R dual e -2^ e ^al = jR -l e 2V> e -$ ) (130) 

$ dual = -il) (131) 

R dual = R _1 exp(-$ - <& dual ). (132) 
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Using (127) to find the transformed dilaton gives 


dual = -$-21n R, 


(133) 


and hence 


R dual = R. (134) 

It can be explicitly checked that equations (27), (34) and (35) are invariant under changes to the dual 
quantities. Equations (27), (28), (34) and (35) provide the integrability conditions for equations (36) and 
(37) which in turn determine the function x- Since the integrability conditions are invariant under the above 
transformation, the equations remain integrable and by substituting the dual quantities into (36) and (37) 
X dual is found to be 


x dual = x _ 2^ + $ + hi R + C 

where C is a constant. 

The dual metric functions (96)-(98) are found as follows 

e c Re®Ai, 

A3, 

A-2■ 


j^aual _ 

a dual _ 

~ 

a dual _ 

— 


(135) 


(136) 

(137) 

(138) 


Since 6(r, r) = 0 in this section, <J?(r, r) satisfies an equation similar to that for i/). To find the corre¬ 
sponding Kasner exponents (cf section 5.1) of the dual model it is necessary to set M = 1 and °tui =° <pi in 
the equations of section (5.1). This reduces to 


dual 

-1 

_ r |+ 0 $i+i( 0 $i) 2 +(°bi) 2 

(139) 

dual 

2 

~ T-2- 0 ^ 

(140) 

dual 

3 

~ T5+V 1 . 

(141) 


Changing to proper time and reading off the Kasner exponents as described in section (5.1) results in 


= 1 

i= 1 

f 2 -1 .. + 

h [1+ °4>i + i(°4>i) 2 + (Vi) 2 ] 2 ‘ 


(142) 

(143) 
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This shows that the behaviors of the original and dual model are very similar. This is expected since 
only Ai really changes and A 2 and A 3 are just exchanged. 

At this point, it should be mentioned that the usual general-relativistic constraint on the sum of the 
Kasner exponents defining the quasi-Kasner behavior is recovered (cf equations (109) and (142)) since we 
are working in the Einstein frame. Assuming a Bianchi I background in the Einstein frame and transforming 
the Kasner solutions from the Einstein to the string frame results in a constraint on the sum of the squares 
of the Kasner exponents being unity. This behavior is characteristic for Kasner-like solutions in the string 
frame [23]. This, in a way, is more illuminating, since it reflects directly the invariance under scale factor 
duality which implies the (discrete) transformation of a Kasner exponent to its negative (pi —> —pi). 

7 Discussion 

We have shown that it is possible to find exact inhomogeneous cosmological solutions of low-energy string 
cosmology. These solutions are cylindrically symmetric and represent cylindrical axionic, dilatonic, and grav¬ 
itational waves propagating inhomogeneously on a flat anisotropic background. When the inhomogeneities 
are of small amplitude these solutions will approach the behavior of small perturbations of isotropic and 
homogeneous anisotropic string models. These solutions also allow us to study the evolution of the universe 
in two physically distinct limits: when the inhomogeneities are larger or smaller than the particle horizon. 
The behavior found has a simple physical interpretation. When inhomogeneities are larger than the horizon 
they evolve quasi-homogeneously but when they enter the horizon there is time for self-interaction to occur 
and the inhomogeneities oscillate like waves. The axion and dilaton fields behave like two fluids in which 
the sound speed equals the speed of light and so shock waves do not form even when the non-linearities are 
of large amplitude. The global structure of our solutions prevents the formation of gravitationally trapped 
regions and so there is no primordial black hole formation. (If the S 3 topology had been chosen, with the 
associated choice R = sinr, then this would have been possible). 

Solutions of varying degree of generality to (3+l)-dimensional string cosmology with dilaton and axion 
in a spacetime of cylindrical symmetry have been discussed. We found that, in general, the axion-dilaton 
system is decoupled from the gravitational background by the cylindrical symmetry. However, the solutions 
of section (3.3) are special in a sense that they describe a universe at large r which contains scalar and 
gravitational waves that are coupled by the wave-like solutions in the axion field. The most general Charach- 
type solutions describe at large values of r an anisotropic universe filled with gravitational and scalar waves 
caused by the dynamics of the axion and dilaton. These two regimes also allow us to find the asymptotic 
behavior of the universe as r -» 0 and r —» 00. There is an initial curvature singularity where the density 
of the dilaton and axion fields is formally infinite (hence we venture outside of the low-energy string theory 
regime assumed here). The early-time behavior resembles the Kasner singularity of general relativity with 
spatially varying indices and is analogous to that observed on scales larger than the horizon at later times. 
The late-time evolution cannot straightforwardly be compared with the present universe because of the 
absence of fermionic fields which provide the standard matter and radiation components of the Big Bang 
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model. The impact of duality upon these solutions is more subtle than in the cosmological models that have 
been examined previously in string theory because of the presence of inhomogeneity. This was discussed in 
detail in section 6 together with the relationships between the results in the Einstein and string frames. 

In summary: we have found exact inhomogeneous and anisotropic cosmological solutions of low-energy 
string theory with non-zero axion and dilaton stresses. These provide a new theoretical laboratory in which 
to explore the ramifications of low-energy string cosmology and to use as a basis for incorporating the effects 
of higher-order corrections. 
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Appendix 

Using the properties of the Bessel and Neumann functions [19] the limits for ip(r, r) (or cj(t, r)) and G(r, r) 
are found. 

The limit nr -C 1 

V’(t, r) ~ a\ (•*/>; r) + a 2 V r ) In r 

where 

00 2 

cci V; r) = Vo + Y] cos [n(r - r n )][ A T n + - In n] 

2 00 

«2 V; r) = Vi + - Y) cos[n(r - r n )] B T n 

7T — , 
n= 1 

The function G(r,r), of eq. (47), is found to approach 

GV;r, r) ~ 7i(V’;r) +7 2 V;r)lnT 

where 

OO 

7i(V’; r) =° V’o +° V’l cos [ n ( r “ r n )] [ A ^n + B ^n In n] 

71= 1 

^ OO 

72 W = -^V’l) 2 +° tpi cos [ n ( r - r n )] B ^n- 
1 n=1 
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The limit nr S> 1 


ip(r, r ) ~ Pi {ip) + P2 {ip) In r + r 2 h(ip ; r, r) 

where 

A WO = Vo 
#2 WO = Vi 

h{ip ; r, r) = V — ) cos[ra(r - r n )][ A ^ n cos(nr - -) + B sin (nr - -)] 
—' \ 7m / 4 4 

n=l 

For G(ip; r, r), the limiting behavior is found to be 

G{ip-r,r) ~ 73 WO +74 WO lnr + 75WO 7 " 


where 

73 (VO =° V’o 

74 WO = ^(Vi ) 1 2 3 4 5 6 7 8 9 

^ 00 

7S V) = ^ICVn) 2 + ( B +n) 2 ] 

n=l 
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